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Abstract 

We investigate the solution of the gap equation for mixed order parameter 
symmetry states as a function of fihing using a two-dimensional tight-binding 
model incorporating second-neighbor hopping for tetragonal and orthorhom- 
bic lattice. The principal (major) component of the order parameter is taken 
to be of the d^2_y2 type. As suggested in several investigations the minor 
component of the order parameter is taken to be of the dxy type. Both the 
permissible mixing angles and tt/2 between the two components are consid- 
ered. As a function of filling pronounced maxima of dx2_y2 order parameter 
is accompanied by minima of the d^y order parameter. At fixed filling, the 
temperature dependence of the two components of the order parameter is also 
studied in all cases. The variation of critical temperature Tc with filling is 
also studied and Tc is found to increase with second-neighbor hopping. 
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I. INTRODUCTION 



Many years after the discovery of the unconventional high- Tc cuprates with a high 
critical temperature Tc the symmetry of its order parameter continues the subject of active 
investigation. From both experimental observations and related theoretical analyses 
]7|^T^ there is a consensus that immediately below the symmetry of the order param- 
eter is of the dx2-y2 type. However, many experiments and theoretical analyses 
p5| , |26|p8| -^ suggest that at lower temperatures the order parameter has a mixed sym- 



metry of the d^-i-y-i + exp(i0)x type, where x represents a minor-component state with a 
distinct symmetry superposed on the major component d^2_y2. From theoretical analysis 
the mixing angle 6 can have the values 0, 7r/2, vr, and 3'7r/2. For mixing angles 7r/2 and 37r/2 
the time-reversal symmetry is broken for these states p5| , p6| . Four possible candidates for 
the minor symmetry state x are the d^y, s, Sx2+y2 and s^y states. From a group theoretical 
point of view these states belong to the same irreducible representation of the orthorhombic 
point group. However, there is still controversy about the nature of the minor component 
and the value of the mixing angle in different high-Tc materials. There have been extensive 
discussions in recent literature about the nature of the minor component and the value of 
the mixing angle. Without repeating those discussions in detail here we present a brief 
summary of these analyses which suggest a mixed-symmetry state and refer the interested 
reader to these references [25,26,28-33,36-50|. 

There are experimental evidences based on the Josephson supercurrent for tunneling be- 
tween a conventional s-wave superconductor (Pb) and twinned or untwinned single crystals 
of YBa2Cu307 (YBCO) that YBCO has a mixed-symmetry state with a minor s-wave com- 
ponent at lower temperatures [15-|T8[1. Kouznetsov et al. performed c-axis Josephson 
tunneling experiments by depositing a conventional superconductor (Pb) across the single 
twin boundary of a YBCO crystal and from a study of the critical current under an applied 
magnetic field they also found evidence of a mixed-symmetry state with a minor s-wave 
component. Similar conclusion was reached by Sridhar et al. ||2^ from the measurement of 
the microwave complex conductivity of high quality YBa2Cu307_5 single crystals at 10 GHz 
using a high-Q Nb cavity. 

From an analysis of angle- resolved photoemission spectroscopy experiment. Ma et al. 
2T| , |22| detected a temperature-dependent gap anisotropy in Bi2Sr2CaCu208+i. which can 
be associated with a mixed-symmetry state with a minor d^y component. More recently, 
Krishana et al. p3[ reported a phase transition to a mixed-symmetry state with a minor 
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dxy component in Bi2Sr2CaCu208 induced by a magnetic field from a study of thermal 
conductivity as a function of temperature and applied field. Later, the experimental study 



of thermal conductivity in Bi2Sr2CaCu208 by Movshovich et al. also leads to the same 
conclusion. 

Different theoretical investigations p5| , |26| provided explanations of the observation by 
Krishana et al. that for weak magnetic field possibly a minor time-reversal symmetry- 
breaking state is induced on a dominating dx2-y2 state in Bi2Sr2CaCu208. This mixed- 
symmetry state is likely to be a minor d^y component on a dx2_y2 state at lower tempera- 
tures. The above investigations in Bi2Sr2CaCu208 suggest a minor d^y component in the 
order parameter. In view of this we perform a theoretical study of this problem using 
a tight-binding model including orthorhombic distortion and second-nearest-neighbor hop- 
ping. Both orthorhombicity and filling are expected to play a crucial role in the evolution 
of these mixed-symmetry states and we study in this paper the properties of the mixed- 
symmetry states of the type dx2_y2 + xdxy at different temperature, filling, orthorhombicity, 
and second- neighbor hopping. A similar study of the mixture of the dx2_y2 and different 
types of s states will be reported elsewhere. 

There is no suitable microscopic theory for high-T^ cuprates and there is controversy 
about a proper description of the normal state and the pairing mechanism for such materials. 
Also, for underdoped systems the existence of a pseudo gap above the superconducting 
critical temperature can not yet be understood within a microscopic theory. However, with 
the increase of doping the superconducting critical temperature increases and the pseudo gap 
is reduced and eventually it disappears at optimal doping. In the absence of a microscopic 
theory, we use a phenomenological two-dimensional tight-binding model with appropriate 
lattice symmetry for studying some of the general features of the mixed-symmetry states 
of overdoped high-Tc superconductors where there is no pseudo gap. This model has been 
used successfully in describing many properties of high-Tc materials [^ , ^ ,^8H35| . 



In high-Tc cuprates, it is generally believed that superconductivity resides mainly in 
the Cu02 planes with nearly tetragonal symmetry. Hence, quite expectedly, the values 
of critical temperature and superconducting gap of high-Tc materials are very sensitive to 
the level of doping which determines the number of electrons available for conduction in the 
two-dimensional Cu02 plane of the cuprates. However, a large anisotropy in the penetration 



depth between the a and h directions has been observed experimentally in YBCO and 
this cannot be understood in a purely tetragonal system. Therefore, in YBCO there is the 
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evidence of a weak orthorhombic distortion. Hence we shall consider orthorhombic distortion 
in the present study of mixed-symmetry states. 

One of the objectives of the present study is to investigate how the superconductivity 
in the uncoupled and mixed-symmetry states changes with doping or filling. The present 
single-band model can acomodate a maximum filling of two electrons (spin up and down) 
per unit cell. In addition to studying the evolution of the mixed-symmetry states with 
respect to filling or doping in presence of second-nearest-neighbor hopping for tetragonal and 
orthorhombic lattice, we also study the temperature dependence of the order parameters for 
a fixed filling. 

In Sec. II we describe the formalism. In Sec. Ill we present our numerical results and 
finally in Sec. IV we give a brief summary. 



We base the present study on the two-dimensional tight binding model which we describe 
below. This model is sufficiently general for considering mixed angular momentum states on 
tetragonal and orthorhombic lattice, employing nearest and second-nearest-neighbor hop- 
ping integrals. In the present tight-binding model the effective interaction V^q for transition 

from a momentum q to k is taken to be separable, and is expanded in terms of some general 
basis functions r/jk, labeled by index i, so that 



The separable nature of the interaction facilitates the solution of the gap equation. The 
orthogonal functions Ty^k are associated with a one-dimensional irreducible representation 
of the point group of square lattice C4„ and are appropriate generalizations of the circular 
harmonics incorporating the proper lattice symmetry. Here Vi is the coupling of effective 
interaction in the specific angular momentum state specified by the index i above. In general 
we have the following 77 's 



II. THEORETICAL FORMULATION 




(1) 



^4q 



V2q 



V5q 



?7lq 



^3q 



cosQx — (3cosqy, 



2 sin g^; sin gj,. 



COS^a; -|- P COS Qy, 

2cosqxCOsqy, 



1, 



s-wave. 



Sxy-wecve, 



dj;2_j,2-wave. 



S2.2_|_y2-wave, 



(2) 
(3) 
(4) 
(5) 
(6) 
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etc. Here (3 = 1 corresponds to square lattice and (3^1 represents orthorhombic distortion. 
The orthogonality property of functions r^'s is given by 

q 

In the present investigation of mixed-symmetry states involving dx2_y2 and dxy, it is ap- 
propriate to consider only two terms in the sum (|ID with rji and above. We shall consider 
singlet Cooper pairing and subsequent condensation in dx2_y2 and d^y waves, denoted by 
indices 1 and 2, respectively, and the mixed-symmetry state formed by these two. Con- 
sequently, the potential couplings Vi and V2 will be taken to be positive corresponding to 
attractive interactions. 

We consider a single tight-binding two-dimensional band with electron dispersion relation 
including second-nearest-neighbor hopping. In this case the quasiparticle dispersion relation 
relating the electronic energy ek and momentum k is taken as 

ek = — 2t(cos kx + P cos ky — 27 cos k^ cos ky) — fi, (8) 

where t and jSt are the nearest-neighbor hopping integrals along the in-plane a and b axes, 
respectively, and 'jt is the second-nearest-neighbor hopping integral. In Eq. (H) fi is the 
chemical potential measured with respect to the Fermi energy and is determined once the 
filling n is specified. The nearest-neighbor hopping parameter t is typically taken to be 
~ 0.1 eV. The parameter j3 destroys the symmetry between the a and b directions in the 
Cu02 planes in this simple model. The potential Vkq above also possesses such a symmetry- 
breaking term. The energy ek is measured with respect to the Fermi surface. Such a 
one-band model with different first-neighbor-hopping parameters in the a and b directions 
is the simplest approximate way of including in the theoretical description the effect of 
orthorhombicity. 

At a finite temperature T, one has the following gap equation 

with i?q = [(cq — /i)^ + lAqp]-*^/^, and k^ the Boltzmann constant. 

It has been observed that the critical temperature Tc is sensitive to the level of doping 
which determines the number of available conduction electrons in the CUO2 plane. In this 
model the chemical potential fi and the the filling n are determined by the number equation 

n = l- V^V^tanh^. (10) 
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The filling n can be related to the the experimental doping S in the three dimensional 
Brillouin zone hj n = 1 — 6. In this work we study the variation of n from to 1 (half 
filling) . 

The order parameter Aq has the following anisotropic form: 

Aq = Air^iq + CA2r/2q, (11) 

where C = exp{i9) is a complex number of unit modulus |Cp = cos^ ^ + sin^ = 1. We 
substitute Eqs. (|I|) and (|ll]) into the gap equation (H) and using the orthogonality property 
(0) obtain the two following coupled equations for Ai and A2: 

In order to proceed further we have to know the nature of C, e. g., whether C is real, 
imaginary or complex. As |Cp = 1, a purely real C implies C = ±1 and in this case no 
further simplification of the coupled Eqs. ( |T2l) and (p^Sf ) is possible and one has 

A V^TA Ai?7iq ± A2r72q , , -Eq ...s 

= ^ym. tanh^, (14) 

±A2 = E V,v,,^^^^^^^±^ tanh (15) 

^ 2Eq 2kBT' ^ ' 

with 

= [(eq - + (Al^lq ± ^2^2^^?^'^ (16) 

The upper sign in Eqs. (0), (P^ , and (|16D corresponds to a mixing angle 6 = Q {C = 1) 
and the lower sign to = vr (C = —1). 

However, Eqs. ([l2|) and ( |l3|) get simplified when C is imaginary, e. g., C = ±i. Conse- 
quently, 

i?q = [(eq-/^)' + A?< + A2<]'/^ (17) 
and Eqs. ( [I^ ) and (|T3|) reduce to 

A.-i:V.^ta.h^, (18) 
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since in this case 

V 2Eq 2kBT ' ^ ^ 

which follows from the definitions of r^iq and ?72q and Eq. (|TTp. However, Eq. (^Dj), which is 
responsible for the simplification for C = ±i {9 = 7r/2 and 37r/2), does not hold for C = ±1 
or for a general complex C. 

Finally, for a general complex C one can separate in a straightforward fashion Eqs. 
(|l^) and (|13D into their real and imaginary parts. In that case Eq. (^) is not valid, and 
this results in four equations for the two unknowns Ai and A2. These four equations are 
consistent only if Ai = or A2 = 0, which means that there could not be any mixing 
between the two components. So mixed-symmetry states are allowed only for mixing angles 
6 = 0, it/2, it, and 37r/2, or for C = ±1, and ±i and we shall consider only these cases in the 
following. Essentially the same conclusion was reached in a different context in Ref. in 
a discussion of mixed-symmetry s-d states on the continuum. In that case mixing was also 



possible only for mixing angles 0,n/2,n and 3tt/2. However, in Ref. |^2[ the discrete sum 
over q was replaced by a continuum integral over angle. We note that the C = ±i cases are 
physically equivalent and lead to the same gap equations ([I8D and (p^Q]). On the other hand, 
C = ±1 cases corresponding to = and vr, respectively, lead to two different sets of gap 
equations (|l^ and ([T5|). 

The ultraviolet momentum-space divergence of the Bardeen-Cooper-Schrieffer equation 
was originally neutralized by a physically- motivated Debye cut off P5| , ^ . This procedure 
had the advantage of reproducing the experimentally observed isotope effect. It can also be 
handled by using the technique of renormalization [|T2|-[l^. Here we introduce a cut off in 
the momentum sums of the gap equation. As there is no pronounced isotope effect in the 
high-Tc cuprates, the present cut off is merely a mathematical one without any reference to 
the phonon-induced Debye cut off. In Eqs. ([T^ ) and ([T3| ) both the interactions Vi and V2 are 
assumed to be energy-independent constants for |eq— /i| < /c^Td and zero for |eq— /i| > kBTu, 
where /c^Td is the present cut off. 



III. NUMERICAL RESULT 



We solve the coupled set of equations (|14D and (|T5D or (|T8D and (|T9D in conjunction 



with the number equation (10). This model gives us the opportunity to study the mixed- 
symmetry state composed of d^i-yi and dxy for different filling and temperature for both 
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tetragonal and orthorhombic lattice. This study will also provide the interesting variation 
of the chemical potential and critical temperature with filling. We solve the above sets 
of equations numerically by the method of iteration and calculate the gaps Ai and A2 at 
various fillings and temperatures. 

We perform calculations on a perfect square lattice with second-nearest-neighbor hopping 
contribution. In addition to the superconductivity in coupled d^2_y2 + d^y and (ij.2_y2 + id^y 
waves we also consider the uncoupled d^^-yi and dxy waves. The effect of orthorhombic 
distortion on our investigation is also studied. The results of our study have interesting 
variation as the second nearest hopping parameter is varied and this is studied in detail 
in the following for mixing angles and 7r/2. Throughout the calculation we consider the 
cut off ksTo = O.lt eV with the parameter t = 0.2586 eV. This corresponds to a cut-off of 
Td = 300 K. The results for the chemical potential /i and order parameters A^.2_y2 and A^y 
presented in this work are all in units of t. 

First, we study the superconductivity in uncoupled dx2_y2 and d^y waves by setting 
A2 = and Ai = 0, respectively, in the equations above at T = on square lattice (/? = !). 
For dx2_y2 wave we use potential parameters Vi = 0.73t and V2 = 0, and for d^y wave we use 
Vi = and V2 = 2. It. The resulting A's with the variation of filling n are plotted in Figs. 
1(a) and 1(b), respectively, for three different values of the second-nearest-neighbor hopping 
parameter 7: 7 = 0, 0.1 and 0.2. In the present study we shall always use Vi = 0.73t on 
square lattice; for 7 = this leads to a critical temperature of 71 K The value of V2 

would be chosen to have a reasonable mixture between the two components. 

From Fig. 1(a) we find that for 7 = 0, Ai = Ax2_y2 has a pronounced maximum at half 
filling and it reduces quickly to zero as n decreases. As second nearest hopping 7 increases, 
the position of this maximum shifts to smaller n. For 7 = 0.2, this maximum is at n ^ 0.8. 
In Fig. 1(b), we see that for 7 = 0, A2 = A^^j, also has a maximum at n = 1; however, this 
maximum is much wider than that in the case of Ai. The position of this maximum also 
shifts towards smaller n as the parameter 7 increases. However, the change in the nature 
of the curves as 7 increases is more pronounced for dx2_y2 in Fig. 1(a) than for d^y in Fig. 
1(b). The variation of A with n in Fig. 1(a) is also quite distinct from that in Fig. 1(b). 

Now we study superconductivity in coupled dx2_y2 + d^y wave at T = on a square 



lattice, which corresponds to mixing angle 0. In this case Eqs. (|lj) and (^) are applicable. 
It is found that with 7 = 0,/? = 1, and Vi = 0.73t, mixing between the two components 
occur for V2 > 1.7t. Here we employ a V2 slightly larger than this critical value to have 
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adequate mixing. The parameters for this model are the following: /3 = 1, Vi = 0.73t, and 
V2 = 2. It. In Fig. 2 we plot Ai = I\^i_yi and A2 = I^xy in this case for different filling 
n and for 7 = 0, 0.05, 0.1, and 0.2. Although the nature of Ai of Fig. 2 is similar to 
that of Ai of Fig. 1(a), the present maxima are narrower. The order parameter A2 in the 
present mixed-symmetry case has a different nature for larger n compared to the uncoupled 
case. Here, in the presence of second-neighbor hopping, as n decreases, A2 first decreases 
rapidly to a minimum and then increases and finally decreases to zero. Moreover, the order 
parameter A2 for this mixed wave has a minimum at the same position of the maximum of 
Ai. 

Next we study superconductivity in coupled d^2_y2 + id^y wave at T = on a square 
lattice, governed by Eqs. (^) and ([ToD , which corresponds to the mixing angle 7r/2 or 37r/2. 
The parameters for this model are the following: /5 = 1, Vi = 0.73t, and V2 = 1.7t. In this 
case to have mixing between the two components for 7 = the critical value of V2 is l.Ot. 
In Fig. 3 we plot Ai = A^.2_y2 and A2 = A^j^ in this case for different filling n and 7 = 0, 
0.05, 0.1, and 0.2. The qualitative nature of Ai of Fig. 3 is similar to that of Figs. 1(a) and 
2. The order parameter A2 for this coupled wave also has a minimum at the same position 
of the maximum of Ai. The minimum in A^.^^ is less pronounced in this case than in the 
case of dx2_y2 + d^y considered in Fig. 2. 

Next we study the behavior of the mixed-symmetry states in the presence of orthorhombic 
distortion (/? = 0.95). First we consider the dr^i^yi + d^y wave. In this case we needed to 



increase V\ to 0.97t in order to have a critical temperature of 71 K for 7 = [Q. The 
parameters for this model are the following: /5 = 0.95, Vi = 0.97t and V2 = 2. It. To have 
mixing in this case the critical value of V2 is 1.5t for 7 = 0. In Fig. 4 we plot Ai and A2 in 
this case for different filling n and second neighbor parameter 7 = 0, 0.05, 0.1, and 0.2. We 
find significant change in the behavior of the order parameters in this case compared to the 
case of square lattice discussed above. In this case A^y could become negative for non-zero 
values of the second-neighbor hopping parameter 7 for values of n close to half filling and the 
phase of the mixing changes with a variation of n in the presence of second-nearest hopping 
7. For a non-zero 7 the mixing angle 6 changes from to vr as n increases. Also, with the 
increase of 7 the zero of A^^^ shifts to a lower value of n, 

Finally, we consider the dx2_y2 + id;j.y wave on orthorhombic lattice with the following 
parameters: (3 = 0.95, Vi = 0.97t and V2 = 1.7t. For 7 = 0, the critical V2 for mixing in 
this case is 1.2t. In Fig. 5 we plot Ai and A2 in this case for different filling n for 7 = 0, 
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0.05, 0.1, and 0.2. Here, the maximum of l^^'^^yi is also accompanied by the minimum 
of IS.xy. The order parameter Ai behaves in a similar way as in Fig. 4, but A2 does not 
become negative. Both in Figs. 4 and 5 we observe a dip in Ai near the maximum. This 
behavior was not found in Figs. 2 and 3 for square lattice. In Figs. 2 — 5 we do not find 
any superconductivity for very low values of filling. 

Now we investigate the temperature dependence of the order parameters in different 
cases. We studied several cases for different values of n, 7 and /3(= 1, 0.95). The qualitative 
nature of the temperature dependence in different cases is quite similar, e. g., the order 
parameters start from a large non-zero value at T = 0. With the increase of temperature 
they first remain essentially constant and then decrease very slowly. With further increase 
of temperature they decrease more rapidly and eventually become zero at a large enough 
temperature. Figure 6 illustrates a typical example of this behavior for /5 = 1, n = 0.8, 
and 7 = 0.2 for the dx2_y2 + d^y and d^^-y^ + idxy cases. The potential parameters in these 
cases are the same as in Figs. 2 and 3. The present temperature dependence of the order 
parameters is quite similar to our earlier studies for n = 1 for dx2_y2 + d^y and dx2_y2 + id^y 

mm- 



waves 



In Fig. 7 we plot the chemical potential /i versus n for the four different models considered 
in Figs. 2 — 5 for the dx2-y2 + dxy and dx2-y2 + id^y cases on tetragonal and orthorhombic 
lattice. The parameters for the model are the same as in Figs. 2 — 5. Surprisingly enough we 
find that fi is essentially the same for the two types of mixing dx2-y2 + dxy and dx2-y2 + idxy, 
independent of the lattice symmetry. However, fi depends on the parameter 7 denoting 
second-neighbor hopping. So in Fig. 7 we have four different curves for four different values 
of second-neighbor hopping. The nature of variation of fi with n changes with 7. The fall 
in /i with n is most rapid in the case of 7 = and slows down as we increase 7. 

In Figs. 8 (a) and 8 (b) we plot Tc versus n for the models of Figs. 2 and 3 for 
three values of the second-neighbor hopping parameter 7: 0, 0.1, and 0.2. We find that 
Tc has a pronounced maxima as a function of n. The maximum value of Tc increases with 
the increase of 7 for both dx2-y2 + dxy and dx2_y2 + idxy waves. A similar increase in the 
superconducting critical temperature with second-nearest-neighbor hopping term was also 
observed in a recent study of temperature and filling dependence of the superconducting 
phase in the Penson-Kolb-Hubbard model ||45[| . 

Finally, we study the interesting ratio 2Aniax/rc as function of filling n in different 
cases. The maximum value Amax of Ak was found numerically in the k plane and the ratio 
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2Aniax/^c was calculated in different function of filling n. As a function of n, 

the quantity Amax has a rapid variation (not shown explicitly here) similar to Tc of Fig. 8. 
Although Tc as well as Amax are found to exhibit rapid variation as function of n, the ratio 
2Aniax/^c is found to be a reasonably smooth function of n. This ratio varies from 4.0 for 
n = 1, to a maximum of about 4.7 as n decreases; the maximum is obtained for n ~ 0.4. 
The average value of 2Aniax/T'c is found to be about 4.3. 



IV. SUMMARY 



In this work we have studied different properties of the mixed drc2_y2+dxy and d^2_y2+idxy 
superconducting states in two-dimensional cuprates using a tight-binding model with appro- 
priate lattice symmetry. Recent experiments and other theoretical investigations suggested 
the formation of such mixed-symmetry states in cuprates on tetragonal and orthorhombic 
lattice. We have included the effect of orthorhombic distortion of the lattice as well as of 
second nearest hopping on the order parameter of these states in the tight-binding model. 
The allowed mixing angles are 0, tt/2,tt, and 37r/2. We studied the variation of order pa- 
rameters A^2_y2 and A^y with filling n for tetragonal and orthorhombic lattices for different 
second-neighbor hopping for pure and mixed-symmetry states. These results are reported 
in Figs. 1 — 5. In the case of mixed-symmetry dx2_y2 +dxy wave on orthorhombic lattice in 
the presence of second- neighbor hopping, the mixing angle changes from to tt as filling n is 
varied. We have also studied the temperature dependence of the order parameter under dif- 
ferent situations. This temperature dependence is quite similar to the n = 1 case studied by 
us previously pl| , p5| . We also studied the dependence of the critical temperature Tc for the 
mixed-symmetry state on filling n under different situations and find that this dependence 
exhibits a maximum and rapid variation. The maximum of Tc increases with the increase of 
second-nearest hopping. The ratio 2Amax/7"c is a slowly varying function of n and varies 
from 4 for n = 1 to a maximum of about 4.7 as n decreases. 

We thank Conselho Nacional de Desenvolvimento Cienti'fico e Tecnologico and Fundagao 
de Amparo a Pesquisa do Estado de Sao Paulo for financial support. 
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Figure Captions: 

1. The order parameters (a) Aj.2_j^2 and (b) A^y at different filling n for second nearest 
hopping 7 = (full line), 0.1 (dashed-dotted fine), and 0.2 (dashed line) for square lattice. 
The potential parameters for (a) are Vi — 0.73i, V2 — 0, and for (b) are Vi — 0,V2 — 2.1t. 

2. The order parameters A^i_y2. (dashed fine), A^y (full line) for different filling n for 
mixed d^^^-yi + d^y wave for second nearest hopping 7 = 0, 0.05, 0.1, and 0.2 on square 
lattice. The parameters for the model are j3 = 1, Vi = 0.73t, V2 = 2. It. 

3. Same as in Fig. 2 for mixed dx2_y2 + id^y wave. The parameters for the model are 
(3=l^Vi = 0.73t, V2 = l.7t. 

4. Same as in Fig. 2 on orthorhombic lattice. The parameters for the model are /3 = 0.95, 
Vi = 0.97t, V2 = 2.1t. 

5. Same as in Fig. 3 on orthorhombic lattice. The parameters for the model are P — 0.95, 
Vi = 0.97t, V2 = 1.7t. 

6. Temperature dependence of the order parameters A^i^y^ and A^y for mixed-symmetry 
dx2-y2 + dxy (dashed line) and d^^-y^ +idxy (full line) waves for (3 — 1, n — 0.8, and 7 = 0.2. 
The larger order parameters represent the dx2_y2 wave. The parameters of the models are 
the same as in Figs. 2 and 3, respectively. 

7. Chemical potential fi versus filling n for the four models of Figs. 2 — 5 for 7 = (full 
line), 0.05 (short dashed line), 0.1 (long dashed line), and 0.2 (dashed-dotted line). Four 
different models lead to identical result. 

8. Critical temperature Tc versus n for the models of (a) Fig. 2 and (b) Fig. 3 for 7 = 
(full fine), 0.1 (short dashed line), and 0.2 (long dashed line). 
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